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Abstract 

The reliability of the mean-field approach to polymer statistical mechanics is investigated 
by comparing results from a recently developed lattice mean-field theory (LMFT) method to 
statistically exact results from two independent numerical Monte Carlo simulations for the 
problems of a polymer chain moving in a spherical cavity and a polymer chain partitioning 
between two confining spheres of different radii, ft is shown that in some cases the agreement 
between the LMFT and the simulation results is excellent, while in others, such as the case of 
strongly fluctuating monomer repulsion fields, the LMFT results agree with the simulations 
only qualitatively. Various approximations of the LMFT method are systematically estimated, 
and the quantitative discrepancy between the two sets of results is explained with the dimin- 
ished accuracy of the saddle-point approximation, implicit in the mean-field method, in the 
case of strongly fluctuating fields. 
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1 Introduction 



Recently, there has been renewed interest ^ in the venerable problem of 

polymer partitioning [|, ^. This phenomenon is of interest both from a theoretical 
as well as a practical perspective. Macromolecular separation techniques, such as 
size exclusion chromatography, gel electrophoresis, filtration, membrane separation, 
etc. all have their basis in polymer partitioning — the property of macromolec- 
ular chains to distribute themselves in a network of random obstacles according to 
their molecular weight/length or electric charge. These phenomena are also of fun- 
damental interest for understanding the properties of macromolecules and complex 
fluids. 

In some recent work p], Q we developed a lattice field theory approach to the 
statistical mechanics of polymers in solution. This approach relied on several approx- 
imations: the treatment was at the mean-field level, which amounts to a saddle-point 
approximation to the system's partition function; also the problem was treated on a 
lattice, converting a continuous three dimensional (3D) problem to a discrete, finite 
matrix problem, and, finally, the polymer chain was assumed to be long enough so 
that it was approximated as a continuous object, ignoring the discrete character of 
the building monomers. 

The goal of this paper is to investigate the range of validity of these approxima- 
tions. In order to do so we need to make a comparison with an exact solution of 
the problem, which can be obtained in the case of neutral polymers by numerical 
Monte Carlo simulations. We are mostly interested in the reliability of the mean- 
field approach and its implicit saddle-point approximation, and the cases in which 
it can be applied with a high level of accuracy. 

In Section 2 of the paper we generalize our previous treatment to the case of a 
soft short-range interaction potential between the monomers in the polymer chain, 
which is needed for the direct comparison with the simulation results. In this section 
we also discuss the other two types of errors, namely the lattice discretization error 
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and the error due to the approximation of the chain as a continuous object — we 
shall call it the "chain continuity error" — and the need for their elimination in order 
to estimate the reliability of the mean-field approach per se. In Section 3 we show 
how the chain continuity error can be eliminated by a transfer matrix representation 
of the partition function, and discuss the effects of this error for the problem of a 
polymer chain in a spherical cavity. Here also we show how the transfer matrix 
approach can be extended to the more general case of a 3D lattice. In Section 
4 we present results for the monomer distribution of a polymer chain in a sphere 
from the transfer matrix approach and compare them with those for the case of a 
continuous chain, and with the exact simulation results, showing a good agreement 
between the three approaches. In Section 5 we describe the procedures used in two 
different numerical simulation methods, and present the results obtained by these 
two independent methods compared with the results obtained by the mean-field 
method using the transfer matrix approach, for the problem of a polymer chain 
confined to move within the region of two connected spheres. The results show that 
for certain range of parameters there is a substantial deviation between the mean- 
field approach and the simulations, which is clearly due to the large fluctuations in 
the monomer repulsion field. We are presently working on a method that would 
allow us to overcome this problem by applying the mean-field approach separately 
in different regions where it is accurate, and "welding" the results for the polymer 
partition function of the separate regions to obtain an accurate description of the 
more complicated physical situation. In Section 6 we outline possible applications 
of the methods and conclude. 
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2 A strategy for studying the accuracy of the mean-field 
approximation in polymer problems 

In order to study systematically the level of accuracy of the mean-field approach to 
polymer statistical mechanics introduced in some recent work of the present authors 
[§, ^, it will be convenient to introduce a slightly generalized Hamiltonian which 
allows for a soft Yukawa type repulsive potential between the monomers of the 
polymer chain [^]. In such a model (absent long-range electrostatic interactions, 
which are not treated in this paper) an essentially exact (up to statistical errors) 
solution of the statistical mechanics is available via standard Monte Carlo simulation 
methods, while, as we shall see below, the mean-field Schrodinger approach of p, |^ 
can be generalized to such a model, allowing a direct comparison and evaluation 
of the error incurred by the saddle-point approximation implicit in the mean-field 
equations. 

The partition sum for a Gaussian polymer chain of M monomers with a repulsive 
Yukawa potential acting between any pair of monomers on the chain and a general 
exclusion potential \4xci(0 (which can be used to exclude monomers from certain 
regions of space) can be written 



^pol 



Here s, s' are integers labelling the monomer location on the chain and the Yukawa 
repulsion between two monomers with coordinate separation r is 

The unconventional choice of ^ for the range of the potential will simplify the 
algebra later on (in particular, the attentive reader will recognize that Vyuk is exactly 
the Green's function of the differential operator 1 — 7A). The partition function 
in (1) involves a positive Boltzmann measure with short range interactions and 
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is therefore susceptible to a standard Monte Carlo simulation. We have found it 
convenient to perform the simulation by updating the monomer locations with a 
heat-bath algorithm for the Gaussian part of the Hamiltonian in (1), followed by 
a Metropolis accept/reject step to include the effects of the Yukawa and exclusion 
potentials Vyuk, Kxci (cf- Section 4 for more details). 

In the mean-field approach to polymer equilibrium statistical mechanics devel- 
oped in our recent work we have generally assumed that the polymer chains 
are sufficiently long that the discrete sums over the monomer index s in (1) may 
be replaced by a continuous integral over a dimensionless variable s, also running 
from to M = total number of monomers in the chain. From the point of view 
of the original model (1) (for which we shall have essentially exact — up to control- 
lable statistical error — results via simulation) this is an approximation, which we 
have termed the "chain continuity approximation" (to distinguish it from the lat- 
tice discretization approximation, for example, in which fields on continuous 3-space 
are replaced by fields defined on a 3D spatial lattice). In this approximation, (1) 
becomes a path integral 



with 

Of) = j ds5{f—x{s)), (4) 

where we have introduced a line distribution j{r) with support along the line tracing 
through the monomers. To establish contact with the models studied in our previous 
work Q , we linearize the dependence on j by introducing an auxiliary field uj via 
a Hubbard-Stratonovich transformation. The path integral (3) then becomes 



^pol 



Dx[s)Duj[r)e > . (5) 



With the range ^/^ set to zero, we recover the simpler model studied in our earlier 
mean-field papers p, |] for the special case of uncharged polymer chains. Here we 
have specialized to the case of a single particle potential V^xci in (§) which is either 



zero or +00, and the sole effect of which is to restrict the fields in the path integral 
to have support in the region where Kxci = 0. The path integral over x{s) 
can be replaced in the standard way by an equivalent Schrodinger problem for the 
imaginary time evolution (from to M) of a particle of mass S/a^ in an imaginary 
potential iXuj{r). Thus 

with 

/ N f r^^, ^ -:r^ jo dsx{s)-i\^ dsu{x(s)) 

Zschr(c^) = / Dx{s)e . (7) 



The mean-field approximation to the system amounts to rerouting the func- 

tional integral over u through a complex saddle point at u = —iuc, where Uc is a 
real function, at which point the evaluation of the Schrodinger amplitude ^schr in 
(0) reduces to a conventional 3D quantum mechanical evolution (in imaginary time) 
of a system subject to the Hamiltonian 

H^-^A + Xu;,{r). (8) 
o 

The hermitian Hamiltonian (^ has a complete set of normalized eigenfunctions 
\l/„(r) with corresponding eigenvalues En, and the evolution amplitude ^schr ap- 
pearing in (0) is just 

^Schr = / dXodXM < Xo|e"*'^^|fAf >= / (ixorfxA,/^^„(xo)^„(xM)e"*^^". (9) 

It has been conventional to assume ground state dominance in evaluating Zschi, but 
this condition is frequently violated in the systems studied below 0, Q . In a recent 
paper, we have presented the general solution |^ which makes no assumption of 
ground state dominance. Defining Zschr = e^p°', and 

An = J «„(r-'), (10) 

we have 

F,o, = \nj:Ale-'''^". (11) 



The saddle point of the functional integral (j^) is then found by minimizing the 
functional 

F = ^l dfucif^il - 7A)^,(r) + Fpoi(a;,). (12) 
Varying this functional as in our previous work @], we obtain: 



where 



Ap(f), (13) 



-ME„ ^-ME„ 



p(r) = -'"^ ^ (14) 

Therefore, 

= (1 - 7A)u;.(f) - p(f) = . (15) 

Equation (|^) implies 

uJcif) = pif) = J Vvukif- r')p{f')dr', (16) 

so that the saddle-point value of the auxiliary field uj{r) can be thought of either as a 
smeared version of the monomer density p or as the total repulsive potential at f^due 
to all monomers. We may therefore eliminate the auxiliary field altogether in (H), 
leading to the following mean-field equation describing the equilibrium properties of 
the polymer chain: 

^V^*„(r-) = A/3(f)^„(f) - E„^„(f). (17) 
o 

This equation is a generalization of Eq. (22) in our previous work 0] for the Yukawa 
type interaction potential between monomers. It is a continuous mean-field nonlinear 
Schrodinger type equation which in general must be solved numerically by putting 
the system on a discrete spatial lattice. This introduces (in addition to the errors 
incurred by making the chain continuous and the saddle-point approximation of (^)) 
a third source of error in comparing the results with the exact solution (via Monte 
Carlo) of the original system (|l]). The determination of the qualitative effects and 
relative size of all three sources of error in some systems of interest is the primary 
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focus of this paper. In particular, it will be important to understand the size of the 
chain continuity error and lattice discretization effects in determining the extent to 
which the mean-field approximation per se fails or succeeds in a specific polymer 
problem. 
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3 Discrete vs continuous chains: a transfer matrix approach 



In this section we estimate the size of the error incurred by replacing a discrete chain 
of M+1 monomers in (|l]) by a continuous imaginary-time Schrodinger evolution as in 
(0). For simplicity, we first study this issue with A = (no intermonomer repulsion) 
for a polymer chain confined to a spherical region of radius R. In this simple case the 
main qualitative effect of this replacement is to force the monomer density to vanish 
exactly at the boundary of the sphere, whereas there is a finite average density of 
monomers at the boundary in the case of the discrete chain |[T^. It turns out that a 
transfer matrix formalism allows a simple analytic description of the discrete chain 
case, which enables us to derive an explicit formula for the size of this "wall effect" . 
Referring to the partition sum ([1|), introduce a transfer matrix T{x, x') as follows: 

T(f,f') = 05(^)e5(x')e-"(^^-^^')', a^T^, (18) 

where Qs{^) is the Heaviside theta-function with support inside the sphere of radius 
R. Viewing the transfer-matrix as a linear operator kernel, the partition function 
(IID can then be written simply as: 

Zpoi = J dxodxMT^ixo,XM)- (19) 

On the other hand, we can define a Hamiltonian ifd for the discrete chain by the 
usual connection, T = e^^"^, so that 

Zpoi = J dxodxM < Xo\e~'^'^^-^\xM > ■ (20) 

Of course, the operator H in (^ and (P) and differ, and it is precisely this 
discrepancy which we have referred to previously as the "chain continuity error". 
Note that the integrals over initial and final monomer locations in Eqs. @ and (|20|) 
(corresponding to free boundary conditions for an open polymer chain) restrict us 
to s-wave (spherically symmetric) eigenf unctions of both H and H^- In the case of 
the Schrodinger operator H, these eigenfunctions are the familiar spherical Bessel 
functions and the corresponding eigenvalues are related in the usual way to the zeroes 
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of these radial functions [as \l/„(r=i?) = in the continuous chain approximation]. 
For the discrete chain operator, we must diagonahze the integral kernel T in (|TB|). 
The eigenvalues e„ of T are related to the eigenvalues En of by = e^^", while 
the eigenfunctions of T and iJd are the same. Restricting ourselves to spherically 
symmetric eigenfunctions, we need to solve the linear integral equation 

g-a|r-r'p^(^/)^^/ _ ^^(^^^ (21) 

r'<R 

which reduces after integration over angles to the one dimensional (ID) integral 
equation 

- r dr' ( e— (^--')' - e-'^(^+'-')') fir') = ef(r), (22) 
a Jo ^ ' 

where fir) = r?/'(r). This equation can easily be solved by discretization and nu- 
merical diagonalization, whence one finds the "wall effect" previously mentioned, 
namely /(-R) 7^ (this effect is also visible in explicit simulations, cf. Section 4). 
We can estimate the size of the effect in the limit <^R (which implies aW' ^ 1 
and the negligibility of the second exponential for r ~ i? in (p2|) ) as follows. To 
leading order in ap/R, the radial transfer matrix in (^)) is 

e-'^ir-r')'^f{r')dr' ^ (l + /(O- (23) 

2 - S 



The ground state corresponds to the eigenvalue — of so in this limit the 
e eigenvalue of the ground state is Sq ~ {l — 4^^)- Next, writing /(r) ~ 

5 + sin (^r-^j , r ~ i?, and substituting this ansatz into (^) at r = i?, we find 

/ e-"('^'-«)' 1^5 + sin 1^"^ ^ dr' ^ ^5. (24) 

The radial integral can be evaluated in the given approximation by linearizing the 
sine function near r = R, sin ~ ^{R — r'), whereupon we find 

— 5~-W-5H . (25) 

71 2 \ a 2aR ^ ^ 

Substituting the value Eq found above for the ground state, we find the size of the 
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wall effect to leading order in l/aE? to be 




(26) 



It should be noted that this estimate is conditioned on /(r) being approximated 
by sin(^7r-^^, that is, /(r) is normalized to be approximately unity at r = R/2. 
In the next section we shall present numerical results for the single sphere system 
confirming this analytic estimate. 

For more general problems in three dimensions, lacking spherical symmetry (such 
as the two sphere situation studied in Section 5) it is still possible to diagonalize 
the transfer matrix, in this case by going to a discrete spatial lattice with lattice 
spacing a. Then the transfer matrix (in this case, truly a finite dimensional matrix!) 
becomes 

Tfirn = e ^> , (27) 

where is the value of the single-particle potential Kxci in (H) at the lattice point 
n (recall that we are still taking A = in (0)). Although in principle a non-sparse 
matrix, we have found it perfectly adequate to cutoff the matrix elements of ( pTf ) for 
\n — m\ > 6, whereupon the remaining relatively sparse matrix can be diagonalized 
by standard Lanczos techniques. The results of calculations employing this transfer 
matrix, and therefore eliminating the chain continuity error in comparisons with the 
Monte Carlo simulations of (|l|), are given in the following two sections. 

If the intermonomer repulsion VVuk is present then the transfer matrix formalism 
can still be employed to estimate the chain continuity error, but only in the context 
of the mean-field approximation discussed in the previous section. After introducing 
the auxiliary field uj and going to the complex saddle point at uj = —iuJc, we can 
reinstate a discrete chain by using the transfer matrix (^) with taken equal to 
the value of Kxci + AcUc at the lattice point n, where Uc is again given by the formula 
(|T6|) and p(r) by ([T^), with slightly different coefficients A^^n defined as 
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4 Comparison of transfer matrix and simulation results for 
a polymer chain in a spherical cavity 



In this section we consider a polymer chain of 100 monomers, restricted to move 
within the volume of a spherical cavity. We have chosen the Kuhn length to be 
unity and a sphere radius equal to ten. First, we discretize the integral equation 
(p2D (neglecting the second exponential) on a fine ID lattice, thus converting it into 
a standard matrix eigenvalue problem, which we solve to a high precision. The 
solution involves only spherically symmetric s-waves. In Fig. 1 we plot the function 
/(r) defined in the previous section by Eq. (p2D, which clearly exhibits the "wall 
effect" discussed earlier — it has a finite value at the edge of the sphere. The value 
f{R) in Fig. 1 is 0.139, which is close to the value calculated from Eq. p6| ) for the 
given set of parameters, 6 ~ 0.145. 

Then we apply the Lanczos procedure described in our earlier work M, B] to 



obtain the eigenvalues and eigenvectors of the transfer matrix (Jig), and therefore 
the monomer density. The matrix is discretized on a cubic lattice of 44 points on each 
side, with lattice spacing a = ap/2, in which the spherical cavity is embedded. The 
resulting radial monomer density is plotted on Fig. 2, where it is compared with the 
exact result from a Monte Carlo heat bath simulation (described in the next section), 
the ID transfer matrix result of Eq. ([2^), and the result from the Schrodinger 
equation approach, Eq. (0), for A = 0. We see that the exact simulation result 
practically coincides with the results obtained via the transfer matrix approach. 
In these cases the monomer density has a finite value at the wall. In the case of 
the Schrodinger approach, where the chain continuity error is present, the radial 
monomer density vanishes at the wall, as expected. 

Since above we have considered a Gaussian chain without excluded volume ef- 
fects, there is no saddle-point error in this calculation; we are left with the lattice 
discretization and chain continuity errors. From a comparison between the ID and 
3D transfer matrix results in Fig. 2, it is clear that in this case the lattice discretiza- 
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Figure 1: The radial function /(r) exhibiting a finite value at the edge of the sphere. 

tion error is practically non-existent, as the results from the fine ID lattice virtually 
coincide with these from the rough 3D lattice. Prom the same figure it is also clear 
that the major error in this case of no excluded volume interactions is the chain 
continuity error, which is also insignificant. 

In Fig. 3 we show an analogous plot of the radial monomer density with all 
parameters the same as in Fig. 2 but with excluded volume interactions present 
with A = 0.0012. It is clear that in this case there is an error due to the saddle- 
point approximation, leading to a deviation between the transfer matrix and the 
simulation result. Here, again, the "wall effect" is present in the simulation and the 
transfer matrix density, while the density resulting from the Schrodinger equation 
approach vanishes at the wall. 
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Figure 2: A comparison between results from the Monte Carlo heat bath simulation (exact), the 
transfer matrix approach, and the Schrodinger equation approach, for the radial monomer density 
P(r) of a Gaussian polymer chain of 100 monomers without excluded volume interactions with a 
Kuhn length unity in a spherical cavity of radius ten. (The error bars of the simulation curve are 
within the size of the symbols used to mark it.) The "wall effect" is less apparent in this figure as 
we are plotting the monomer density (roughly speaking, the square of /(r) in Fig. 1). 



14 




Figure 3: A comparison between results from the Monte Carlo heat bath simulation (exact), the 
transfer matrix approach, and the Schrodinger equation approach, for the radial monomer density 
P{r) of a polymer chain of 100 monomers with excluded volume interactions with A ~ 0.0012. The 
Kuhn length is unity and the spherical cavity has a radius of ten. (The error bars of the simulation 
curve are within the size of the symbols used to mark it.) 
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5 Comparison between lattice field theory and simulation 
results for a polymer chain partitioned between two con- 
nected spherical cavities of different radii 

The simulation of the system (|l|) by Monte Carlo techniques is readily accomplished 
given the absence of long-range forces. We have chosen to perform the simula- 
tion by generating new polymer configurations using a heat-bath algorithm on the 
Gaussian part of the Boltzmann measure e " , rejecting the obtained 

configuration if it is forbidden by the exclusion potential Vexci (here assumed to be 
either zero or oo), and then correcting for the remaining Yukawa potential with a 
Metropolis accept/reject step. Specifically, the algorithm employed is as follows: 

1. First, a monomer location i is chosen at random (i.e. an integer between and 
M — 1) and a new location Xi generated according to the weight 



e oc e-3(^^-*) /"p, (29) 

where di = |(fi-i + Xj+i). 

2. If the new location Xi is outside the allowed region (as determined by the 
exclusion potential) it is rejected and a new Xi found by the preceding algorithm. 

3. The Yukawa part of the Hamiltonian, Sy^k = f J2s,s' VYuk( used to 
accept /reject the new configuration (Metropolis update), with Yukawa Hamil- 
tonian S'yuk as follows: 

(a) if ^Yuk < 'S'Yuk, the new configuration is accepted without further ado, 

(b) if S'Yuk > "S'Yuk, the new configuration is accepted with probability e~'^Yuk+'5'Yuk_ 
If it is rejected, the monomer is left unmoved. 

In order to obtain confidence in our results, we have performed a second dy- 
namic Monte Carlo simulation based on the kink-jump technique developed by 



Baumgartner and Binder |[TT|. In this method the polymer chain is treated as a 
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Figure 4: A polymer chain confined to move within two connected spheres. 

"pearl-necklace" in which every two consecutive monomers are connected by a rigid 
rod of length /. Monomers are modeled as hard spheres, whose radius 7*772 can De 
adjusted to describe the strength of the excluded volume interaction. Initially the 
chain is placed in the allowed region with the positions of its monomers chosen ran- 
domly. Then the chain dynamics is evolved by the kink-jump technique: at each 
time step a monomer is picked, say the nth one, and then it is rotated around the 
axis connecting the [n — l)th and the (n + l)th monomers at a random angle. If 
an end monomer is chosen, it is moved to a new random position keeping the rod 
length between it and the monomer to which it is connected fixed. The move is 
accepted if the monomers do not overlap and stay within the allowed region. 

We apply these two simulations methods to a polymer chain confined to move 
within two connected spheres of different radii. In relative units, we choose the big 
sphere to have a radius -Ri equal to unity, the small one has a radius R2 = 0.8, the 
Kuhn length is chosen to be Op = 0.2, and the spheres are embedded in each other 
with the distance between their centers fixed at 1.7. This system is schematically 
shown in Fig. 4. 
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We have calculated the polymer partition coefficient K, defined as the ratio of 
the average number of monomers in the big and small spheres respectively, K = 
(Ml) / (M2), where the total number of monomers is M = (Mi) + (M2), for varying 
M. First, in Fig. 5 we show the simulation results for the case of a non- interacting 
chain. These results are compared with lattice field calculations using the transfer 
matrix and the Schrodinger equation approaches. As in the case described in the 
previous section, we see that the agreement between all the different approaches is 
very good for the case of zero interaction between the monomers, which confirms 
once again that the chain continuity and lattice discretization errors are small. In 
this figure the result from the Schrodinger equation approach appears to be closer 
to the statistically exact simulation result than the transfer matrix result, probably 
due to a coincidental cancellation of errors. Here, and in the next two figures, we 
observe that in the case of longer chains the transfer matrix approach produces a 
partition coefficient which is systematically smaller than the one resulting from the 
simulations. 

In Fig. 6 we show results for In K from the simulations for the interacting chain 
with Ai = 0.0012 and 7 = for system (||), which is matched to a monomer radius 
'^m,! = 0.025 in the kink-jump technique. These results are also compared to the 
results obtained via the transfer matrix lattice field calculation for the same param- 
eters used in the Monte Carlo heat bath/metropolis calculation. In this calculation 
the two-sphere system is modeled on a discrete cubic lattice of 44 points on each 
side with lattice spacing a = 0.1 as described in [|, In our previous work |^ 
we have shown that the size of the aperture connecting the two spheres plays an 
important role in determining the value of the partition coefficient K. However, 
when the system shown in Fig. 4 is modeled on a lattice, lattice artifacts prevent us 
from obtaining the exact same aperture size, as in the continuous case. Therefore, 
the aperture size on the lattice is adjusted to be as close as possible to the one in 
the continuous case shown in Fig. 4. 



18 



5.0 



0.0 



- 


1 ' 1 ' 1 ' 1 ' — 


1 ^ 

T - 








o 






o ^ ^ 

— 1 — ^ Q 








J^^.y ° heat bath 
^7^^^ ■ kink— jump 

♦ — ♦ transfer matrix 








V ^-v Schr. equation 






t 


I . I . I . I . 







) I 1 I 1 I 1 I 1 I 1 I 1 

0.0 100.0 200.0 300.0 400.0 500.0 600.0 

M 



Figure 5: IniiT vs M for the case of a Gaussian chain with no excluded volume interactions. The 
simulation results are seen to agree with the lattice field calculations. 

In our previous work P] we have shown that the excluded volume parameter A 
can be connected to the monomer radius via the second virial coefficient, leading to 
the expectation that 

A const X r^. (30) 

From the parameters of the calculation whose results are presented in Fig. 6 
we calculate const and use it to calculate a second excluded volume parameter, 
A2 = 0.0033, corresponding to the monomer radius = 0.035, and then we per- 
form simulations for the two new sets of parameters by the heat bath/metropolis 
method and the kink-jump technique, and present the results in Fig. 7. The good 
agreement between this second set of simulations confirms formula (PD|). We have 
also performed the mean-field transfer matrix calculations for this new A, and have 
plotted the results for comparison in the same figure. 

From Figs. 6 and 7 we observe that Ini^ increases almost linearly with M for 
small M, then goes through a turnover, and at large M decreases to a limiting value 
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Figure 6: \nK vs M for A = 0.0012, which is matched to rm = 0.025. The mean-field transfer 
matrix results for the same A are seen to disagree substantially with the statistically exact results, 
although in both curves there is a maximum with a turnover. 



20 



3.0 



2.0 



1.0 



o heat bath/metropolis 
■ kink— jump 
♦ mean— field 



0.0 ' ' ' ' ' ' ' ' ' ' ' ' ' 

0.0 100.0 200.0 300.0 400.0 500.0 600.0 

M 

Figure 7: Ini^ vs M for A = 0.0032928, corresponding to = 0.035 in the kink-jump simulation 
technique. The mean-field transfer matrix results for the same A are seen to disagree with the 
simulation results, although both curves possess the same generic features. 
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determined by the ratio of the volumes of the two spheres. Prom the simulations we 
observe that, as is to be expected, the turnover occurs when the end-to-end distance 
of the polymer chain becomes comparable to the diameter of the larger sphere. Both 
Figs. 6 and 7 show that although the mean-field results show the same generic be- 
havior, they deviate substantially from the exact simulation results. Since we have 
eliminated the chain continuity error by using the transfer matrix approach, and we 
have already seen that the lattice discretization error in these calculations is rela- 
tively small, we conclude that the quantitative discrepancy between the statistically 
exact simulation results and the mean-field results is due mainly to the saddle-point 
approximation implicit in the mean-field approach. Also, it should be noted that 
in these cases the mean-field results from the Schrodinger equation approach (not 
shown in the figures) are very close to the transfer matrix results, leading to the 
conclusion that the chain continuity error is small in the case of excluded volume 
interactions too. 

The results show that the error is larger in the cases of shorter chains, while in the 
cases of longer chains it is relatively small. This can be explained with the help of 
the following observation made in the course of the simulations: the shorter chains 
tend to be localized entirely in one sphere or the other (thus achieving a higher 
conformational entropy), and only make infrequent jumps between the two spheres. 
This creates a highly fluctuating monomer repulsion fleld and renders the saddle- 
point approximation less reliable. In the case of a long chain, monomer repulsion 
competes with the free energy gain due to the higher conformational entropy, and 
the chain tends to stretch between the two spheres, occupying both of them simul- 
taneously for most of the time. In that case the monomer repulsion field exhibits 
smaller fiuctuations and the saddle-point approximation is expected to be better. 
To illustrate this point, in Fig. 8 we show a plot of the number of monomers in 
the bigger sphere. Mi, as a function of the number of simulation sweeps, for two 
different cases, M = 500 and M = 100, for the parameters corresponding to Fig. 7. 
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Figure 8: The number of monomers in the bigger sphere, Mi, as a function of the number of simu- 
lation sweeps t ioT M = 100 (left plot) and M = 500 (right plot), for the parameters corresponding 
to Fig. 7. 

It is seen that the relative fluctuations in Mi are much larger for the smaller M. 

It is also important to point out that what we actually measure in the simulations 
is the average number of monomers in the big sphere, (Mi), and the partition 
coefficient K computed from it is very sensitive to small errors in (Mi), as (Mi) and 
(M2) are constrained to add to M, hence a small error in (Mi) will automatically 
be reflected in (M2) with an opposite sign, thus amplifying the error in the ratio 
K, especially when {Mi) is much bigger than (M2), as in the calculations presented 
here. The actual error between the mean-field calculations and the simulations is, 
in fact, much smaller than one would infer from Figs. 6 and 7. To illustrate this 
point, in Fig. 9 we plot the results for (Mi) and (M2) as a function of M from the 
calculations corresponding to the In K plot in Fig. 7. 



23 



o simulation 
■ mean-field 



0.0 100.0 200.0 300.0 400.0 600.0 
U 



o simulation 
■ mean-fiekl 



200.0 400.0 
U 



Figure 9: The average number of monomers in the bigger and smaller spheres, (Mi) (left plot) 
and (M2) (right plot), respectively, as a function of the total number of monomers, M, for the 
parameters corresponding to Fig. 7. 
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6 Conclusions 



The main focus of this work has been investigation of the vahdity and rehabihty of 
mean-field methods in polymer statistical mechanics by systematic elimination and 
estimation of all errors inherent to the lattice field theory approach developed in our 
earlier work ^, apart from the saddle-point approximation implicit in the mean- 
field method. To be able to evaluate these errors we have compared the lattice field 
theory results to the statistically exact results from two independent Monte Carlo 
simulations, namely, a heat bath/metropolis computation and a computation based 
on the kink-jump simulation technique. In the case of no excluded volume interac- 
tions, computation of the radial monomer density of a polymer chain constrained to 
move within a spherical cavity, and of the partition coefficient of a chain restricted 
to move within two connected spheres, shows that "chain continuity" and "lattice 
discretization errors" are relatively small. On the other hand, including short-range 
pair-wise interactions in the Hamiltonian and subsequent mean-field treatment of 
these via the saddle-point approximation leads to a larger quantitative error in the 
partition coefficient for problems of certain geometry, as is shown by a comparison 
with the exact simulation results. 

This error is larger for shorter chains, and diminishes for long chains. It has been 
attributed to the strongly fluctuating monomer repulsion fleld in the case of short 
chains, as revealed by the simulations: the short chains preferentially localize entirely 
in one sphere or the other (thus attaining a higher conformational entropy), making 
infrequent jumps between the two spheres, which leads to a strongly fluctuating 
repulsion fleld felt by each monomer in the chain, and diminishes the validity of 
the saddle-point approximation, an approximation which assumes smoothly varying 
flelds. In the case of long chains, the large number of monomers leads to a very 
strong repulsion, overwhelming the gain in conformational entropy which would 
have been achieved had the whole chain been localized in one of the two spheres, 
hence the chain stretches into both spheres in its most likely conflguration, creating 
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a smoother interaction field felt by the individual monomers. Thus, in this case the 
saddle-point approximation is expected to be good, and indeed, the mean-field and 
the simulation results show good agreement. 

The calculations presented in this paper suggest that the mean-field approach 
works very well for problems of polymers moving easily in an open region, but is 
not very accurate in situations where the polymers are forced to move between sev- 
eral such regions connected only by narrow conduits leading to large fluctuations in 
the monomer density field. Fortunately, in certain cases it turns out to be possi- 
ble to "weld" the mean-field results for the polymer partition function in separate 
open regions into an accurate description of the more complicated situation, so that 
the saddle-point approximation is only involved where it is accurate. In this way 
we have been able to reproduce the simulation results of Figs. 6 and 7 for In K. A 
detailed description of this technique will be presented in a forthcoming publication. 
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